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Abstract. We consider binary dispatching problem originating from ob- 
ject oriented programming. We want to preprocess a hierarchy of classes 
and collection of methods so that given a function call in the run-time 
we are able to retrieve the most specialized implementation which can be 
invoked with the actual types of the arguments. This problem has been 
thoroughly studied for the case of mono dispatching [7I4| . where the 
methods take just one argument, resulting in (expected) O(loglogm) 
query time after just linear preprocessing. For the binary dispatching, 
where the methods take exactly two arguments, logarithmic query time 
is possible [5] , even if the structure is allowed to take linear space [1 J . Un- 
fortunately, constructing such structure requires as much as (expected) 
<9(m(loglogm) 2 ) time [Tl9] . 

Using a different idea we are able to construct in (deterministic) linear 
time and space a structure allowing dispatching binary methods in the 
same logarithmic time. Then we show how to improve the query time 
to just 0( 1 J°fJ^ m ), which is easily seen to be optimal as a consequence 
of some already known lower bounds if we want to keep the size of the 
resulting structure close to linear. 

Key-words: method dispatching, persistent data structures, rectangle 
geometry 

1 Introduction 

The motivation for the method dispatching comes from object-oriented pro- 
gramming languages, where we have a hierarchy of classes with uniquely defined 
parents. We also have a collection of m functions accepting a constant number 
of arguments, where each argument must have a specified class as an ancestor 
in the hierarchy. Then given a function call, we should (efficiently) determine 
the most specific implementation based on the actual types of the arguments. 
This problem was first considered in the mono dispatching version, where each 
method takes just one argument. It is known that in such case the input can 
be preprocessed in linear time and space so that each query can be answered in 
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Fig. 1. Unique lowest bridge for u and v and an ambiguous situation. 

O (log log m) time [7], and that we can update the structure in C(m e ) time while 
retaining the same bounds on the query time [5]. 

The more general multi-method version of the problem was considered by 
Ferragina et al. [5] , whose methods imply that for the special case of the binary 
dispatching (where all function are binary) we can achieve O (log log m) query 
time after a 0(m 1+£ ) preprocessing or C(log m) query after a 0(m log m) prepro- 
cessing. Then Eppstein and Muthukrishnan j3j improved the former by showing 
that, for example, 0(1) query time is possible after a 0{m 1+t ) preprocessing. Fi- 
nally Alstrup et al. pQ decreased the preprocessing space to 0(m) while retaining 
logarithmic query time. Unfortunately, their preprocessing time was not linear 
but (expected) C(m(loglogm) 2 ). Another structure with the same bounds was 
given by Poon and Kwok [9]. 

In this paper we first give in Section [3] a simpler yet more effective solution 
with (deterministic) linear time and space preprocessing and the same logarith- 
mic query time. Our solution uses a slightly different approach than the previous 
methods, and because of this difference we are then able to decrease the query 
time in Section 4 to just O ( i^°fj^ m ) while retaining the linear time and space 
preprocessing. This complexity is easily seen to be optimal for structures of size 
O{m\og c m) as a consequence of some already known lower bounds, which we 
briefly review in Section [5j 

While even the first logarithmic query time solution needs the word RAM 
model, the same is true for the previously known linear space solutions, hence it 
should not be seen as a drawback. 

2 Preliminaries 

We work with the following formulation of the problem: we are given a tree T on 
n vertices and a collection of m bridges, which are simply pairs of vertices (it, v). 
We say that a bridge (u, v) is lower than (u' , v') if u' is a descendant of u and v' 
is a descendant of v. We want to preprocess the input so that given two vertices 
u' and v' we can detect the lowest bridge («, v) such that v! is a descendant of u 
and v' is a descendant of v, see Figure [lj If there is no such unique lowest bridge, 
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Fig. 2. Rectangles on the same stack and their version tree. 



we need to signal ambiguity. We aim to develop a 0(n + m) time preprocessing 
which allows O(logm) time queries. 

We work in the standard RAM model of computation with logarithmic word 
size. The following result is known in such model. 

Lemma 1 (atomic heaps |6j). It is possible to maintain a collection of sets 
S(i) C {1, 2, . . . , n} so that inserting, removing and finding successor in each of 
those sets work in constant time (amortized for insert and remove, worst case 
for find) as long as \S(i)\ < log c n for all i, assuming O(n) time and space 
preprocessing, where c is any (but fixed) constant. 



3 New algorithm 

First observe that the above problem reduces in a natural way (by computing 
the pre- and post-order numbers) to retrieving the smallest rectangle containing 
a given point on a n X n grid (or detecting there is no such unique smallest 
rectangle). From now on we will work with this simple geometric formulation. 
Note that the x and y projections of any two rectangles are either disjoint or 
contained in each other (we call such collection of rectangles valid) and we can 
normalize the coordinates so that n < 2m. 

We sweep the grid from left to right while maintaining a structure describing 
currently intersected rectangles. The structure is simply a full binary tree on 
n leaves corresponding to different y coordinates. To process an interval [2/1,2/2] 
with j/i < 2/2 , we locate the lowest common ancestor v of the leaves corresponding 
to 2/1 and 2/2 and call it responsible for [2/1,2/2]- Each inner vertex stores a stack 
containing all intervals it is currently responsible for. To insert a new interval 
we push it onto its responsible vertex stack. To remove an interval, locate the 
responsible vertex and observe that (because the collection is valid) the interval 
we want to remove is its top element, and we can simply pop it. 

Fix an inner vertex and consider all I rectangles it was responsible for, see 
Figure [2] Note that the intersection of all their y projections is nonempty, and 
hence any two of those projections are contained in each other. By a simple 
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linear time transformation we can assume that all their start and end points 
are different, and their sorted list is x\ < x 2 < ■ ■ ■ < xi (of course we cannot 
assume that sorting a single list can be performed in linear time, but note that 
we can sort lists of all inner vertices at once, and because their elements are 
small integers we can apply counting sort). We define the version tree as follows: 
the parent of a rectangle [xi,^] x [2/1,2/2] is the rectangle [x'^x^] x [2/1,2/2] such 
that [xi,X2] C [x'^x^] and [x'^x^] is the smallest possible. Because any two x 
projections are either disjoint or contained in each other, and all Xi are different, 
this is a valid definition. We may assume that the result is indeed a tree (not a 
forest) by adding one artificial rectangle. Each vertex of this tree is labeled with 
an integer denoting the height 7/2 — 2/1 of the corresponding rectangle. Consider 
the sorted list of all different x coordinates. For each pair of consecutive integers 
Xi < Xi + i on this list we would like to find the vertex of the version tree such 
that its ancestors are exactly the elements of the stack at time t £ (xi^Xi+i) 
(we call it the tail at time t). This can be precomputed in a straightforward way 
during the sweep. 

Consider a query concerning a point (x, j/). First we locate the leaf v corre- 
sponding to x. Any rectangle containing (x, y) belongs to the stack of one of its 
ancestors at time x. More specifically, it must be an ancestor of the tail at time x 
of one of those log to stacks. Hence we should start with locating all those log to 
tails efficiently. 

Lemma 2. Given a time t we can retrieve its tail at every ancestor of the leaf 
corresponding to v in total O (log to) time after a linear time and space prepro- 
cessing. 

Proof. A straightforward application of the fractional cascading technique of 
Chazelle [2]. Recall that this technique allows linear time and space preprocess- 
ing of a constant-degree graph with (sorted) lists of elements associated to the 
vertices so that given a path we can perform binary search for the same value in 
all lists corresponding to its vertices in time O (log to + p), where m is the total 
size of all lists and p is the length of the path. In our case p = log to and the 
claimed running time follows. □ 

Each version tree will be carefully preprocessed as to implement two opera- 
tions. Let path(i>) be the set of all ancestor weights of a given vertex v. The first 
operation is very simple: given v we would like to find the vertex corresponding 
to maxpath(i;). This can be trivially preprocessed in linear time and space. The 
second operation is more involved: given v and x we would like to find the vertex 
corresponding to the successor of x in path(v). Before we show how to implement 
it efficiently, we formulate two auxiliary lemmas. 

Lemma 3. A collection of sets of points S(i) on anx n grid such that \S(i)\ < 
log 2 n for all i can be preprocessed in 0(n + J2i 1^1) time so that given i and 
(x,y) we can retrieve the point corresponding to mm{y' : (x',y') 6 S(i) A x' < 
x A y' > y} in constant time. 
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Fig. 3. Reducing the original query to queries in smaller structures. 




Proof. The idea is to recursively build a collection of smaller structures allowing 
performing the same operation but on subsets of the original set of point. 

Assume we have a set of m < log 2 n points (xi,yi). By Lemma [I] we can 
sort the points according to their y coordinates in 0(m) time. Then we split 
the original set into blocks of size roughly ^/m by choosing yfm evenly spaced 
elements (x iv ^, in this sorted sequence and call the i-th block Bi. We build 

a smaller structure for each Bi . Additionally, let x\ be the smallest x coordinate 
in the i-th block and create a new set of y/m points of the form (x^, i) which we 
call representatives. We build a smaller structure for this new set. If m is very 
small, say m < \/\ogn, we switch to a different method: we can first normalize 
the coordinates of the points so that they are at most \/Iog m and encode the 
whole set in a single machine word going row-by-row. 

Observe that the total size of all structures built for a single S(i) is just 
0(S(i)). Furthermore, they allow us to answer a single query in constant time as 
follows. First locate the block y belongs to and query the corresponding smaller 
structure. If this smaller structure contains a point (x',y') with x' < x and 
y' > Vi we are done. Otherwise we use the smaller structure built for the rep- 
resentatives to locate the lowermost block containing such point and query its 
corresponding smaller structure, see Figure [3] If the size of S(i) is small and we 
have the whole set encoded in a single machine word, we can answer a query by 
first masking out all bits corresponding to points with too big x or too small 
y coordinates and then finding the lowest bit set to 1. The total running time 
is constant because we will inspect just a constant number of structures for a 
single query. □ 

The next lemma will be used to preprocess each version tree. The main tool 
in its proof is the heavy path decomposition, which is defined as follows: each 
vertex chooses an edge leading to a child with the largest size. Removing all 
non-chosen edges leaves us with a collection of paths. We define the path tree by 
creating one vertex for each such path, and choosing the parent of a path p by 
looking up the parent of its highest vertex in the original tree and retrieving the 
corresponding path. It is easy to see that the depth of a path tree is just logn. 
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We say that a path p is above a vertex v if the path v belongs to is a descendant 
of p in the path tree. 

Lemma 4. A node weighted weighted tree on n vertices with the weights from 
{1, 2, . . . , n} can be preprocessed in linear time and space so that given v and x 
we can find the vertex corresponding to the successor of x in path(w) in O(logn) 
time. 

Proof. Consider a single path. By preprocessing all paths at once we can con- 
struct a sorted list of all weights weight(t>i) < weight^) < ■ ■ • < weight (i^) on 
this path. We split each list by choosing 1q ^ n evenly spaced weights weight(w Q , log 2 „) 
and call the corresponding vertices important. The a-th group contains vertices 
u Qlog 2 n _ A for A = 0, 1, . . . , log 2 n — 1. We choose the highest vertex from each 
such group and calling it the representative. Note that the total number of both 
important vertices and representatives is at most lo ™ 2 - ■ For each group we con- 
struct the corresponding small set: 

(depth(u a i og 2 n _ 4 ),weight(v alo g2„_ 4 )) for all A = 0,1,.. . ,log 2 n- 1 

and apply the preprocessing described in Lemma [3j 

We would like to implement the following operations efficiently: 

1. given v and x, for each path above v find the successor of x among the 
weights of its important vertices, 

2. given v and x, for each path above v find the successor of x among the 
weights of all representatives which are ancestors of v. 

First lets see how such information allows us quick retrieval of the successor of x 
in path(u). Assuming that we know the successor of x among all weights of the 
important vertices on p, we query the structure constructed for the corresponding 
group. If it contains at least one vertex above v, we are clearly done. Otherwise 
we know that the successor belongs to a higher group than x. Assuming that 
we know the successor of x among the weights of all representatives which are 
above v, we query the structure constructed for his group. 

We use almost the same to implement both operations. Lets start with the 
former. For each path p we build a binary search tree containing all important 
vertices located on all paths corresponding to the ancestors of p (including p 
itself) in the path tree. The vertices are sorted according to their weight. By 
using any persistent balanced search trees we can build the structures in total 
linear time and space. Furthermore, the total number of new nodes created as 
a result of all inserts will be just 0{^-^). For the sake of concreteness, assume 
that we use trees in which the original elements are stored only in the leaves. To 
facilitate efficient query processing, at each vertex v we store an additional helper 
structure mapping a path depth to the smallest weight stored at the subtree of v 
and originating from an important vertex with such path depth. This structure 
consists of an array of size log n and one word with the i-th bit set if and only if 
the i-th entry in the array is defined. The helper structures are of just C(logn) 
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size, hence we can afford to build one for each new node in total linear time 
and space. Now given a query concerning a vertex v, we locate its path and 
the corresponding binary search tree. Then wc find the successor of x in this 
tree with a single O(logn) time transversal. We claim that the helper structures 
stored at all right brothers of the visited vertices give us enough information to 
locate the successors of x among the chosen weights of all paths above v. This 
is fairly obvious if we consider a single such path. The tricky part is to extract 
the information from all of them in O(logn) time. We go through the vertices 
in a bottom-up order. In the very beginning we do not have the successor on 
any path. We iteratively consider the right brother of the next visited vertex: its 
helper structure gives us the successor on every path for which the corresponding 
entry is defined and which is yet unknown. By storing the currently unknown 
set in a single word, we can compute this intersection in constant time, and then 
extract the successors in constant time per path. 

To implement the latter, we use almost the same method. The only exception 
is that now we build a binary search tree for each representative instead of each 
path. □ 

Now we are ready to prove the main lemma in this section. 

Lemma 5. A set of m rectangles on a n x n grid can be preprocessed in linear 
time and space so that given a point (x,y) we can find the rectangle [xi,X2] x 
[yi 12/2] containing (x,y) with the smallest height y2 — yi in O(logm) time. 

Proof. First apply Lemma [2] to locate the tail at time x in every ancestor of 
the leaf corresponding to x. Then select v to be the lowest of those ancestors 
such that the maximum on the corresponding tail-to-root path is sufficiently 
big to contain y in the corresponding interval. There are just log to ancestors 
and extracting each maximum requires constant time. Now we claim that the 
smallest height rectangle can be found by looking at the version tree of v only. 
Assume otherwise, i.e., there is v' such that v' is an ancestor of v and the smallest 
rectangle belongs to the version tree of v' . But then the y interval at v' properly 
contains the the y interval at v, and hence the interval of the rectangle at v is 
smaller. 

To finish the proof, we apply Lemma [4] to each version tree. Then given 
the tail at v we first compute the smallest possible height of a rectangle stored 
in this tree which guarantees containing y inside. This can be performed in 
O(logm) time if we store all y intervals sorted according to their lengths. Then 
we compute the successor of this smallest possible height on the tail-to-root path. 
It corresponds to the smallest height rectangle. □ 

Theorem 1. Bridge color problem can be solved in O(logm) time after a linear 
time and space preprocessing. 

Proof. First we transform the input so that no two bridges share an endpoint. 
This can be ensured by increasing the number of vertices to at most n' = n + 2m 
by repeating the following procedure: given a group of bridges (u, Vi), (u, v%) ,. . ., 
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(u,Vf~) with the same endpoint u, sort them so that depth(ui) < depth(«i + i) for 
alii = 1, 2, . . . , k — 1. Then replace u by a path U\ — ?> U2 —>...—> and create 
fc bridges (u,, «j) for i = 1, 2, . . . , fc. We interpret the resulting set of bridges as 
a collection of rectangles on a n' x n' grid. 

By Lemma [5] we can find the rectangle [xi,^] x [3/1 > 2/2] containing (x,y) 
with the smallest height 1/2 — 2/1 in O(logm) time. By swapping all x and y 
coordinates, we can also find the rectangle with the smallest width X2—X\. Note 
that because no two bridges share an endpoint, those two rectangles are uniquely 
defined. If they are different, we signal ambiguity. Otherwise this rectangle is the 
unique lowest bridge. □ 



4 Decreasing the query time 

In order to reduce the query time we will carefully modify the algorithm from 
Section [3) We use a combination of standard tricks (increasing the degree of 
the tree to log 6 n) and a few extension of the observations from the previous 
section. The high level idea stays the same: we sweep the plane from left to 
right maintaining a structure describing currently intersected rectangles. The 
structure is a full tree of degree d = log 6 n for a sufficiently small e (to be fixed 
later) on n leaves corresponding to different y coordinates . Each vertex stores 
three structures: 

1. left stack storing prefixes of the corresponding segment, 

2. right stack storing suffixes of the corresponding segment, 

3. block stack storing fragments consisting of a number of whole segments cor- 
responding to a contiguous range of its children. 

To process an interval [2/1,2/2] we locate the lowest common ancestor u of 
the leaves corresponding to 2/1 and 2/2 and split the interval into three parts: a 
suffix of a segment corresponding to some child of v\ , a number of full segments 
corresponding to t^+i, . . . ,Vj-i (which we call the middle part) and finally a 
prefix of a segment corresponding to some child Vj , where vt, i>2, • • • , f <j are the 
children of u. All stacks will be implemented using technique similar to the one 
from Lemma|4j with the block stack requiring one additional detail. Nevertheless, 
it also uses the idea of storing a list of all different x coordinates where a new item 
appears or disappears. For each pair of consecutive Xi < Xi+i we store a pointer 
to the corresponding version of the structure, and apply fractional cascading to 
quickly locate the most up-to-date version given a query (x,y). Now we need a 
slightly stronger version, though. 

Lemma 6. Given a time t we can retrieve the pointers to all current structures 
at all ancestors of the leaf correspondinq to x in total 0( , losm ) time after a 
linear time and space preprocessing. 

Proof. We apply the fast fractional cascading of Shi and JaJa [TU]. It allows 
a linear time and space preprocessing of a log 6 m-degree (where e < |) tree 
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with (sorted) lists of elements associated to the vertices so that given a path of 
length p we can find the successors of a given value in all lists corresponding to 



The implementation of a block stack consists of two parts. For each i < j we 
store a stack of intervals which the middle part corresponds to v i7 ...,Vj. Note 
that if we store a pointer for each i < j, the space usage might be too large. 
Fortunately, only a linear number of those pointers will be non-null during the 
whole execution, hence for each vertex we can store a dictionary mapping 
to a pointer. The dictionary can be implemented efficiently using LemmaJT] The 
second part of the implementation is a single word encoding information which 
stacks are currently nonempty. Observe that the stack of intervals stored for each 
i < j is indeed a stack: we either push a new interval or pop the one which is on 
the top. Hence its implementation will be the same as in the case of left and right 
stacks. The only difference is that given k we need to quickly find i < j such that 
the corresponding stack is nonempty, i < k < j and j — i is smallest. This can be 
easily done in constant time using the single word encoding all nonempty stacks. 
To finish the implementation we need to show how to implement all stacks. As 
in the previous section, we will store their version trees, and a stack is actually a 
pointer to the current vertex in such tree. Each version tree will be preprocessed 
using a stronger version of Lemma [4] For that we first need to extend Lemma [3| 

Lemma 7. A collection of sets of points S(i) on a n x n x n grid such that 
\S(i)\ < log 3 n for alii can be preprocessed in 0(n + \Si\) time so that given 
i and (x,y,z) we can retrieve the point corresponding to min{z' : (x',y',z') G 
S(i) A x' < x A y' < y A z' > z} in constant time. 

Proof. We extend the method used in the proof of Lemma [3] We partition S(i) 
into yz, xz and xy blocks of size roughly m 3 / 4 by choosing to 1 / 4 evenly spaced 
elements in the sequence of points sorted according to the first, second, and 
third coordinate, respectively. For each such block we build a smaller structure. 
Additionally, we create a new set of representatives by taking all original points 
and replacing their coordinates with the number of the corresponding yz, xz, 
and xy block. Note that the size of this new set is at most m 3 / 4 . We build a 
smaller structure for this set of representatives. If m < log 1 / 3 n, we normalize 
the coordinates and encode the whole set in a single machine word instead. 

To answer a query concerning (x, y, z) we first use the smaller structures built 
for yz, xz, and xy blocks. This gives us the answer if its coordinates are close to 
x, y or z. Otherwise we use the smaller structure built for the representatives, 
which gives us the z coordinate of the answer. Having this coordinate, we use 
the smaller structure built for the corresponding xy block. Overall, the running 
time is constant, and the total size of all structures is 0(m), as the depth of the 
recursion is constant. □ 

To extend Lemma |4j we need a relaxation of the heavy path decomposition, 
which we call a thin fragments decomposition. First we choose all edges connect- 
ing vertices whose subtrees are of size at least j^z—- After removing this top 




its vertices in time 0( 
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Fig. 4. Decomposition of a tree into thin fragments. 



fragment of the tree we get a collection of smaller trees. For each of them we do 
the same, namely we choose all edges connecting vertices whose subtrees are of 
size which is a logarithmic fraction of the whole size of the current tree, remove 
the resulting top fragments, and repeat the whole procedure until we get single 
vertices. Each removed fragment is a tree on at most log n leaves (but possibly 
many more inner vertices, hence the name), see Figure|4] We define the fragment 
tree with vertices corresponding to fragments and edges defined in the natural 
way, namely we make one vertex a child of another if the root of the former frag- 
ment is a child of a vertex belonging to the latter fragment. It is easy to see that 
its depth is at most n ■ Fragment depth is the depth of the corresponding 

vertex in the fragment tree. 

Lemma 8. A node weighted weighted tree on n vertices with the weights from 
{1,2, ... ,n} can be preprocessed in linear time and space so that given v and x we 
can find the vertex corresponding to the successor of x in path(i>) in 0( lo 1 °j 5 O g n ) 
time. 

Proof. For each fragment containing of I vertices we construct a sorted list of 
all weights. We choose lo Jj - evenly spaced weights out of them and call the 

corresponding vertices important (notice log 3 n instead of log 2 n) . For each such 
important vertex v a log 3 n we construct a small set of all pairs 

(-pre(t; Q , log 3 n _ /i ),post(u alog 3 n _ A ),weight(u alog 3 n )) for A = 0, 1,.. . ,log 3 n 

and apply the preprocessing described in Lemma [7J Additionally, from each 
group we select vertices which are minimal in the relation of being an ancestor 
and call them the representatives. Note that a single group cannot have more 
than log n such representatives, therefore the total number of representatives is 

We will show how to preprocess the tree so that given v and x, for each 
fragment above v we can locate the successor of x among the weights of all 
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important vertices and the weights of all representatives which are above v. 
We use exactly the same method as in Lemma |4j namely for each fragment / 
we build a binary search tree containing all important vertices on all fragments 
corresponding to ancestors of / in the fragment tree (including /), with a smaller 
structure stored at each vertex. Similarly, for each representative we build a 
binary search tree containing all representatives above. Now the size of the helper 
structure is just O ( lo '° fj^ n ) , though, hence we can find the successor for each 

fragment in O ( lo ^g „ ) time after just linear preprocessing. Given the successor, 
we use the helper structure storing the whole group of the corresponding vertex 
to retrieve the answer in constant time per fragment above v. □ 

This gives us all the ingredients necessary to speed up the method from the 
previous section. 

Theorem 2. Bridqe color problem can be solved in O( , lo f m ) time after a 

a r v log log m ' J 

linear time and space preprocessing. 

5 Lower bound 

It turns out that O ( ig°f^ m ) query time is optimal if we want to keep the size 
of the structure close to linear. This follows from the results of [8], where a 
lower bound for the following 2D stabbing problem is shown: preprocess a given 
collection of m 2D rectilinear rectangles so that we can quickly retrieve (any) 
rectangle containing a given point. It turns out that if we want to keep the size of 
the structure C(mlog c m), the best possible query time is f?( lo ^" m ), even if we 
allow randomization. Furthermore, the lower bound is shown through a reduction 
from the reachability oracle problem, and hence the x and y projections of any 
two rectangles in the collection are either disjoint or contained in each other. 
A closer inspection of the proof shows that the queries can be assumed to be 
chosen so that there is at most one rectangle containing the point. It follows that 
we can encode the instance as a binary method dispatching problem, and the 
claimed lower bound follows. 

6 Conclusions 

We presented a time-and-space optimal solution for the binary method dispatch- 
ing problem. Two questions remain: 

1. is it possible to make the structure dynamic? 

2. can we achieve linear preprocessing and logarithmic query in the pointer 
machine model? Note that the construction of Alstrup et al. [T) makes a 
heavy use of word RAM specific structures, such as the van Emde Boas 
trees. 
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